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ABSTRACT 

We revisit the Galactic chemical evolution of oxygen, addressing the systematic errors 
inherent in classical determinations of the oxygen abundance that arise from the use 
of one dimensional hydrostatic (ID) model atmospheres and from the assumption of 
local thermodynamic equilibrium (LTE). We perform detailed 3D non-LTE radiative 
transfer calculations for atomic oxygen lines across a grid of 3D hydrodynamic STAG¬ 
GER model atmospheres for dwarfs and subgiants. We apply our grid of predicted line 
strengths of the [OI] 630 nm and OI 777 nm lines using accurate stellar parameters 
from the literature. We infer a steep decay in [O/Fe] for [Fe/H] > —1.0, a plateau 
[O/Fe] ss 0.5 down to [Fe/H] ss —2.5 and an increasing trend for [Fe/H] < —2.5. 
Our 3D non-LTE calculations yield overall concordant results from the two oxygen 
abundance diagnostics. 

Key words: radiative transfer — line: formation — stars: atmospheres — stars: 
abundances — Galaxy: abundances — methods: numerical 


1 INTRODUCTION 

Oxygen abundances and abundance ratios are crucial for un¬ 
derstanding the evolution of stars and galaxies. It is there¬ 
fore disconcerting that the [O/Fe] 1 vs [Fe/H] relationship in 
the metal-poor halo is still in dispute, in spite of tremendous 
effort over a number of years (for an in-depth review, see for 
example Stasinska et al. 2012). 

The situation is complicated by the availability of differ¬ 
ent oxygen abundance diagnostics, that often yield discor¬ 
dant results. Abundances inferred from the low-excitation 
forbidden [OI] 630 nm line tend to show a plateau at 
[O/Fe] R3 0.5 below [Fe/H] « —1 (Nissen et al. 2002; Cayrel 
et al. 2004; Garcia Perez et al. 2006), as do those inferred 
from OH infrared (IR) vibration rotation lines (Balachan- 
dran et al. 2001; Melendez & Barbuy 2002). Analyses of the 
high-excitation permitted OI 777 nm lines (Carretta et al. 
2000; Nissen et al. 2002; Garcia Perez et al. 2006) and of 
the OH ultraviolet (UV) excitation lines (Israelian et al. 
1998, 2001; Boesgaard et al. 1999; Gonzalez Hernandez et al. 
2010) tend to predict increasing oxygen abundances toward 
lower metallicities, with [O/Fe] « 0.7 at [Fe/H] « —2.0 

* E-mail: anish.amarsi@anu.edu.au 

1 The logarithmic abundance of an arbitrary element A is defined 
with respect to hydrogen: log ca = log (//- + 12. The abundance 
ratio of elements A and B is given by: [A/B] = (log e \ — log e®) — 
(log eg — log eg), where © denotes the solar value. 


and [O/Fe] ~ 1.0 at [Fe/H] « —3.0. Notably, however, 
Garcia Perez et al. (2006) inferred a plateau at [O/Fe] « 0.5 
from the OH UV lines in metal-poor subgiants between 
[Fe/H] = —1 and —3. It has been suggested (e.g. Israelian 
et al. 2001) that these different trends could reflect in part 
different rates of circumstellar mixing in dwarfs, subgiants 
and giants, because at low metallicities the different diagnos¬ 
tics are typically used in different types of stars. Spite et al. 
(2005) compared abundances in giants which have and have 
not undergone internal mixing. They found no significant 
systematic difference in [O/Fe], indicating that the problem 
must in fact arise from elsewhere. 

The discrepancies between the [OI] 630 nm line and the 
OI 777 nm lines likely have at least two points of origin. 
Firstly, classical abundance analyses are prone to system¬ 
atic modelling errors. These errors arise from the use of one 
dimensional (ID) hydrostatic model atmospheres and from 
the assumption of local thermodynamic equilibrium (LTE). 
Although the [OI] 630 nm line forms in LTE (Kiselman & 
Nordlund 1995), the effects of inhomogeneities and three- 
dimensional structures in late-type stellar atmospheres on 
the [OI] 630 nm line in the metal-poor regime are significant 
(Nissen et al. 2002; Collet et al. 2007). The OI 777nm lines 
are a classic case study of non-LTE line formation (Sedlmayr 
1974; Kiselman 1993; Takeda 2003; Fabbian et al. 2009). 
However, full 3D non-LTE calculations have not been been 
carried out beyond the Sun (Kiselman & Nordlund 1995; 
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Asplund et al. 2004; Pereira et al. 2009; Prakapavicius et al. 
2013; Steffen et al. 2015). 

Secondly, systematic errors can be present in the 
adopted stellar parameters. In particular, offsets in the 
adopted effective temperature (T e a) would reveal them¬ 
selves as offsets in the [O/Fe] vs [Fe/H] trend, as the 
[OI] 630 nm line and the OI 777 nm lines have opposite sen¬ 
sitivities to T e ff. If T^ff is underestimated for a given star, 
the [OI] 630 nm line will underestimate the atmospheric oxy¬ 
gen abundance, while the OI 777 nm lines will overestimate 
it. The accurate calibrations of the infrared flux method 
(IRFM) by Casagrande et al. (2010) have shown previous 
Teff estimates dwarfs and subgiants with [Fe/H] < —2.5 to 
typically be ~ 200 K too low. 

Offsets in T e ff of this magnitude are not sufficient to 
explain by themselves the [O/Fe] vs [Fe/H] problem at the 
lowest metallicities (Fulbright & Johnson 2003). Likewise, 
extremely large 3D non-LTE abundance corrections would 
be required in the OI 777 nm lines to explain the prob¬ 
lem, because the non-LTE abundance corrections in these 
lines (Fabbian et al. 2009; Dobrovolskas et al. 2014) go in 
the same direction as the 3D abundance corrections in the 
[OI] 630 nm line (Nissen et al. 2002; Collet et al. 2007). Thus, 
errors in the stellar parameters must be addressed simulta¬ 
neously with ID LTE errors to achieve reliable results from 
the oxygen abundance diagnostics. 

In this letter we aim to pin down the true nature of the 
Galactic chemical evolution of oxygen. To that end, we have 
performed for the first time detailed 3D non-LTE radiative 
transfer calculations for atomic oxygen lines across a grid of 
3D hydrodynamic STAGGER model atmospheres (Collet et al. 
2011; Magic et al. 2013) (§2). We illustrate the results in 
the form of 3D non-LTE vs ID LTE abundance corrections 
(§3). Using our 3D non-LTE grid, accurate stellar parame¬ 
ters from the literature, and a large collection of published 
equivalent widths of the [OI] 630 nm and OI 777 nm lines, 
we revisit the [O/Fe] vs [Fe/H] relationship (§4). 


2 CALCULATIONS 

2.1 Model atmospheres 

3D hydrodynamic model atmospheres were taken from the 
STAGGER-grid (Collet et al. 2011; Magic et al. 2013). These 
were trimmed by removing the very optically thick lay¬ 
ers of the original simulations and resampled: the mesh- 
size was thus reduced from 240 2 x 230 gridpoints down to 
120 2 x 110, the last dimension representing the vertical. One¬ 
dimensional (ID) hydrostatic model atmospheres were com¬ 
puted using an equation of state and a treatment of opac¬ 
ities and radiative transfer fully consistent with the STAG¬ 
GER models. 

We illustrate the extent of the grid in effective tem¬ 
perature (Teff) and surface gravity (logg) space in Fig. 1. 
Grid nodes are spaced in intervals of Aloggr = 0.5 and 
AT ef f « 500 K. The grid spans —3.0 < [Fe/H] < 0.0, in 
equal intervals of A [Fe/H] = 1.0. 


OI 777nm [OI] 630nm 
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Figure 1 . Abundance corrections across the STAGGER grid of 
3D stellar atmospheres. Values beyond the periphery of the grid 
were obtained by linear extrapolation and should be treated 
with caution. Abundance corrections are shown across T e ff-log g 
planes. Top: [Fe/H] = 0.0; log eo = 8.7. Bottom: [Fe/H] — 3.0; 
logeo = 6.2. The ID microturbulence is £ = l.Okms -1 . 


2.2 Spectral line formation 

Multi-level statistical equilibrium radiative transfer calcula¬ 
tions were performed using the domain-decomposed MPI- 
parallelised code multi3d (Leenaarts & Carlsson 2009). For 
a given oxygen abundance and 3D model, the average line 
profile was calculated from four temporal snapshots. The 
monochromatic average non-LTE to LTE ratio was then 
applied to the average line profile from at least 50 tempo¬ 
ral snapshots calculated in LTE using SCATE (Hayek et al. 
2011 ). 

The 3D and the ID spectral line formation calculations 
were carried out in the same way to avoid incurring sys¬ 
tematic errors. In ID an additional parameter, the micro¬ 
turbulence £, is required to describe the line broadening by 
convective velocity fields. The ID grid was calculated for 
several values of £; in this letter we adopt a typical value of 
£ = l.Okms -1 to illustrate the abundance corrections (§3). 
We emphasise that in 3D no microturbulent velocities are 
necessary (Asplund et al. 2000). 

2.3 Model atom 

A 23-level model oxygen atom was used, based on the mod¬ 
els used by Carlsson & Judge (1993), Kiselman (1993), and 
Fabbian et al. (2009). The model was updated with energies 
and oscillator strengths from NIST 2 (Kramida et al. 2012), 
natural line broadening coefficients from VALD3 3 (Piskunov 
et al. 1995; Kupka et al. 2000), and collisional line broaden¬ 
ing coefficients from Barklem et al. (1998). 

The statistical equilibrium is sensitive to the adopted 

2 http://www.nist.gov/pml/data/asd.cfm 

3 http: //vald. astro.univie.ac.at/~vald3/php/vald.php 
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collisional rate coefficients. Results from R-matrix calcula¬ 
tions were used for the rate of excitation via collisions with 
electrons (Barklcm 2007). Following Lambert (1993), a cor¬ 
rection to the original formula of Drawin (1968, 1969) was 
adopted for the rate of excitation via collisions with neutral 
hydrogen. All levels were coupled in this manner, including 
including those with radiatively weak and forbidden tran¬ 
sitions, in which case Drawin’s formula was used with an 
effective oscillator strength / m i n = 10 _li . An enhancement 
factor Sh to the Drawin formula was calibrated by study¬ 
ing the centre-to-limb variation in the Sun, following Pereira 
et al. (2009) using observational data from the Swedish So¬ 
lar Telescope. We have obtained Sh ~ 1 and have inferred 
a solar oxygen abundance logeo ~ 8.68. 

The dominant systematic errors in the non-LTE calcula¬ 
tions likely arise from the errors in the rate of excitation via 
neutral hydrogen collisions. As discussed in Barklem et al. 
(2011), the approach of calibrating the Drawin formula us¬ 
ing a single parameter Sh cannot correct any errors in the 
relative rates, and can mask other deficiencies in the mod¬ 
els. One can get a feeling for the sensitivity of the results 
on the neutral hydrogen collisions by comparing the results 
with those obtained when Sh is reduced to nought, thereby 
neglecting them altogether. When doing so, the derived non- 
LTE abundance corrections for the OI 777 nm can become 
up to 0.05 dex more negative. 


2.4 Background opacities 

While background continuum opacities were calculated in¬ 
ternally by the two spectral line formation codes, back¬ 
ground line opacities had to be included separately. These 
were pre-calculated assuming LTE and the chemical compo¬ 
sition of the atmosphere being modelled (Collet et al. 2005; 
Bergemann et al. 2012). 

It can be tempting to simplify the non-LTE analysis 
by neglecting the background lines that do not overlap with 
the lines being studied in detail. This is dangerous; oxygen 
proves to be a case in point. We have found the Lyman-a line 
to play a critical role in controlling the statistical equilibrium 
of oxygen in the photosphere. Without the Lyman-a line, 
extremely large abundance corrections would erroneously be 
inferred from the OI 777 nm lines in metal-poor stars. We 
discuss this briefly in §3, and in more detail in a forthcoming 
paper. 

The [OI] 630 nm line is heavily blended with a Nil line 
(Lambert 1978; Allende Prieto et al. 2001). This Nil line was 
included in LTE in both 3D and ID, with a nickel abundance 
that was scaled relative to the solar value logem = 6.20 
(Scott et al. 2015). 


3 ABUNDANCE CORRECTIONS 

Abundance corrections A log eo = log e ^ D ’ NLTE _ log ^d+te 
were calculated by demanding that the equivalent width of 
the line flux obtained in LTE from the ID model with a 
given oxygen abundance log £q D ’ lte matches that inferred in 
non-LTE from the 3D model with some oxygen abundance 
loge Q ’ . Equivalent widths were determined by direct 

integration across the line profile. We illustrate corrections 



3D NLTE 

Figure 2 . Abundance corrections A logeo = l°g £ o ~ 

logeQ D,LTE for the [OI] 630 nm line ( red circles) and for the 
OI 777 nm lines ( black diamonds) in a dwarf (solid-, T e g = 5771 K, 
log g = 4.44), a subgiant (dash; T e gf = 5000 K, log g = 3.0), and a 
turn-off star (dot-dash; T e g = 6500 K, I rig g = 4.0), with ID mi¬ 
croturbulence £ = l.Okms -1 and varying chemical composition. 
The ID LTE oxygen abundance is 8.7+ [Fe/H] + [a/Fe], where the 
enhancement is [a/Fe] = 0.5 for [Fe/H] < —1.0 and [a/Fe] = 0.0 
for [Fe/H] = 0.0. 


for the two oxygen abundance diagnostics across different 
models in Fig. 1 and Fig. 2. 

The 3D corrections for the [OI] 630 nm line are close to 
zero in the solar metallicity case, but become increasingly 
more negative at lower [Fe/H], in line with previous findings 
(Nissen et al. 2002; Collet et al. 2007). Differences in the 
predicted mean temperature stratification between the hy- 
drodynamical 3D models and the hydrostatic ID models are 
more severe in this regime; ID models typically overestimate 
the atmospheric temperature in the upper layers where this 
line forms (Asplund et al. 1999). 

In contrast, the 3D non-LTE vs ID LTE abundance 
corrections for the OI 777 nm lines are most severe in the 
solar metallicity case, reaching up to A log eo ~ —0.5 (Fig. 1 
and Fig. 2). Photon losses in the line itself drive the non- 
LTE effects (Asplund et al. 2004); inhomogeneities in the 
3D model atmospheres exacerbate these effects. 

Fabbian et al. (2009) obtained extremely large non- 
LTE effects at low [Fe/H] in their ID non-LTE study of 
the OI 777 nm lines. These authors found that, toward 
lower metallicities, photon pumping (Asplund 2005) in the 
OI 130 nm lines would drive very large overpopulations with 
respect to LTE in the upper level of these lines. Then, 
collisional coupling would propagate this overpopulation 
onto the lower level of the OI 777 nm lines. This would in¬ 
crease the OI 777 nm line opacity, strengthening the lines 
and causing large negative abundance corrections, reach¬ 
ing A logeo ~ —0.5 in metal-poor turn-off stars (assum¬ 
ing Sh = 1). However, the analysis of Fabbian et al. (2009) 
did not include background opacity from the Lyman-a line, 
which is significant because the line is very strong and the 
OI 130 nm lines sit on its red wing. We have found that the 
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Figure 3. Final results based on equivalent widths from the lit¬ 
erature (Carretta et al. 2000; Fulbright & Johnson 2003; Bensby 
et al. 2004; Akerman et al. 2004; Garcia Perez et al. 2006; Fabbian 
et al. 2009; Nissen et al. 2002, 2014; Bertran de Lis et al. 2015) 
and stellar parameters from Nissen et al. (2014) and Casagrande 
et al. (2010, 2011). Top: [O/Fe] vs [Fe/H] inferred from the 
[OI] 630 nm line (red circles) and from the OI 777 nm lines (black 
diamonds). Lines of best fit to the data in the domains —2.5 < 
[Fe/H] < —1.0 and —1.0 < [Fe/H] < 0.0 are overdrawn. The fits 
were obtained by minimising \ 2 and assigning each star equal 
weight; where abundances from both diagnostics were available, 
the mean value was used. Bottom: differences in [O/Fe] for indi¬ 
vidual stars inferred from the [OI] 630 nm and OI 777 nm lines, 
for the cases where both sets of equivalent widths were avail¬ 
able. The median value is 0.02 dex and the standard deviation is 
0.14 dex. 


Lyman-a line provides an efficient alternative destruction 
route for UV photons at low [Fe/H] that completely stifles 
the photon pumping effect in the OI 130 nm lines. Thus, in 
the absence of an alternative non-LTE mechanism, we are 
left with small departures from LTE in the OI 777 nm lines 
at low [Fe/H]. 


4 THE GALACTIC CHEMICAL EVOLUTION 
OF OXYGEN 

The equivalent widths of the [OI] 630 nm and 
OI 777 nm lines were taken from a number of studies 
of dwarfs and subgiants based on high signal-to-noise 
observations (Carretta et al. 2000; Fulbright & Johnson 
2003; Bensby et al. 2004; Akerman et al. 2004; Garcia Perez 
et al. 2006; Fabbian et al. 2009; Nissen et al. 2002, 2014; 
Bertran de Lis et al. 2015). Where more than one value 
was available for a given star and line, the unweighted 
mean was adopted. For consistency, equivalent widths 
across the theoretical 3D non-LTE grid were computed by 
fitting Gaussian functions to the line fluxes and integrating 
analytically. Stellar parameters were taken from several 
recent studies (Casagrande et al. 2010, 2011; Nissen et al. 
2014). These T e fr estimates were either derived using or 
critically compared to the accurate IRFM calibrations of 


Casagrande et al. (2010); where more than one set of stellar 
parameters was available for a given star, the newest set 
was adopted. 

Hitherto, studies have typically found discrepant results 
from the two abundance diagnostics at low [Fe/H] (§1). 
In Fig. 3 we compare the [O/Fe] ratios inferred from our 
analyses. We have found the [OI] 630 nm line and the 
OI 777 nm lines to give similar [O/Fe] vs [Fe/H] trends 
down to [Fe/H] « —2.2, the lowest metallicity in which the 
[OI] 630 nm line is detectable in halo subgiants and turn-off 
stars. Furthermore, the abundances inferred from these di¬ 
agnostics in the atmospheres of the same stars are consistent 
to within a standard deviation of 0.14 dex. 

It bears repeating that there are two factors in our anal¬ 
ysis that are absent in most previous studies, that conspire 
to give concordant results between the two abundance di¬ 
agnostics. First, we have accounted for 3D non-LTE effects 
in the OI 777 nm lines: these are of decreasing importance 
towards lower [Fe/H] , but even then remain significant. Sec¬ 
ond, we have used new and accurate stellar parameters: the 
more reliable IRFM calibrations give T e ff estimates that are 
significantly larger than those typically used in the past. 

The [O/Fe] vs [Fe/H] relationship in Fig. 3 reflects the 
evolution with time of oxygen and iron yields. Oxygen is syn¬ 
thesised almost entirely in massive stars (M > 8 Mq), its 
most abundant isotope 16 O being the endpoint of Helium 
burning (Woosley et al. 2002; Clayton 2003; Meyer et al. 
2008). Since iron is synthesised in type II supernovae ex¬ 
plosions (Woosley et al. 2002), the plateau at [O/Fe] « 0.5 
between —2.2 < [Fe/H] < —1.0 indicates that massive stars 
at this epoch eject an approximately constant ratio of oxy¬ 
gen and iron upon their deaths, in agreement with Galactic 
chemical evolution models (Frangois et al. 2004; Kobayashi 
et al. 2006). The steep linear decay seen above [Fe/H] > —1.0 
is usually interpreted as a sign of the long-lived type la 
supernovae becoming significant, increasing the rate of en¬ 
richment of iron into the cosmos (Tinsley 1979; McWilliam 
1997). 

At [Fe/H] < —2.5 there is a slight upturn in [O/Fe], 
but it is less pronounced than found from the ID LTE anal¬ 
yses of OH UV lines by Israelian et al. (1998, 2001) and 
Boesgaard et al. (1999). The upturn could indicate a shift 
in the mass distribution of stars at earlier epochs; stars that 
are more massive and more metal-poor are expected to yield 
larger oxygen-to-iron ratios (Kobayashi et al. 2006). We cau¬ 
tion however that inferences in this region are less reliable 
because they are based on a very small sample of stars and 
on analyses of the OI 777 nm lines alone. A larger sample 
of halo turn-off stars with accurate stellar parameters and 
very high signal-to-noise spectra will be needed to confirm 
this result. 
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